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SUMMARY 

The present methods available for the determination of the 
magnitudes of the vibrational relaxation times of molecules arc 
discussed and the recently developed theory of Schwartz, Slaweky, 
and Herzfeld is used to compute the variation with temperature of 
the collisional vibrational excitation probabilities of mixtures of 
N2 and C»2» and of and HF. The results are used to estimate 
the extent of vibrational temperature lag in the hydrogen -fluorine 
rocket motor. 
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I. INTRODUCTION 

The internal energy, E, of a gas may be written, approxi- 
mately, as the sum of energies of translation, rotation, vibration, 
and electronic motion. That is to say, we write that 

E - f * f„tCT r ) ♦ E,;.CT V ) . e„a e ), (l) 

where T.T^ ,T V , and T e represent, respectively, the absolute temp- 
eratures associated with the translational, rotational, vibrational, 
and electronic energy states. Under some conditions these energy 
states will be in equilibrium such that T*T|.-T^ -T e . If, in the flow 
of a gas, temperature changes occur at rates sufficiently fast to ex- 
ceed the time required for these internal energy states to readjust 
themselves, the gas will depart from its equilibrium partition of 
energy. 

The term "vibrational relaxation time" has. come into use as 
a measure of the time required to establish equilibrium with respect 
to the internal vibrations of the molecules. The other internal states 
adjust themselves so rapidly that we can usually neglect their relaxa- 
tion times. The vibrational relaxation time, T , is defined by the 
relation 

-- T-[E vib (T)-F^(r v )]. ( 2 ) 

Working fluids such as steam, air, and exhaust gases have 
considerable vibrational heat capacity at high temperatures. If these 
gases have relaxation times comparable with or shorter than the in- 
tervals during which temperature changes occur in the gas, looses 
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( 1 ) 

must be expected. For turbine working fluids, Kantrowitz* ' esti- 
mates that in some cases the losses at high temperatures due to this 
adjustment lag can. be comparable with the losses due to skin fric- 
tion. Unfortunately, there exist3 a scarcity of data on the magni- 
tudes of the relaxation times of most gases and, without these values 
being known at various temperature levels, their effects cannot be 
analyzed. 

Various experimental methods have been devised to measure 
relaxation times of gases. Of these, only the two discussed here 
have been used with any degree of success. A third method involv- 
ing a shock tube has been proposed, but has as yet not actually been 



used. 
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H. EXPERIMENTAL METHODS 

A. THE SOUND VELOCITY METHOD* 2 * 

The sound velocity method can be adapted for the determina- 
tion of the relaxation times of gase3 up to temperatures of the order 
of 230Q°K. The two essentials are (a) the production of sound waves 
of a definite frequency and (b) the accurate measurement of temper- 
ature. 

13 ) 

The apparatus, originated by Pierce' ' , consists essentially 
of a carbon tube which can be heated to a high temperature and in 
which a stationary train of sound waves can be set up from a quartz 
crystal vibrating piezo- electrically. A piston traverses the tube. 

The position of the piston is read on an accurately graduated steel 
rule. This piston functions as a reflector to set up a standing wave 
system in the tube. The amplitude of vibration of the crystal changes 
as the reflector is moved along the tube. Resonance is detected by 
observing the magnitude of the plate current of the oscillator which 
maintains the quartz crystal in vibration. The half wave-length is 
determined by recording the different positions of the reflector at 
which resonance occurs. The temperature of the gas in the tube is 
measured by sighting an optical pyrometer of the "disappearing fila- 
ment" type on to the face of the movable reflector. An electric fur- 
nace surrounds the carbon tube so that the temperature may be ad- 
justed to higher levels. Figure 1 shows the general layout of the ap- 
paratus. 

The velocity of sound, V, in the gas is computed from the 
measured half-wave length in the tube and from the known frequency 
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of vibration of the crystal. The ratio of specific heats is given by 
the formula, 

v-V*M/RT, (3) 

v/here M is the molecular wdight, R io the gas constant, and T is 
the absolute temperature. The correction necessary to the above 
formula for departure from the perfect ga3 laws' is negligible, since 
at the temperatures involved most gases are essentially perfect 
gases. From the observed data, an apparent ra.tio of specific heats, 
which we may call y ', is computed from Eq. (3). Also from the 
relationship, 

C v /R - , <4) 

an apparent specific heat at constant volume, C.^, can be obtained. 

The crystal used for the production of sound waves, if cut 
properly, can be made to vibrate in either the flexural mode or the 
longitudinal mode by suitable electronic circuitry. The two fre- 
quencies produced may vary by as much as 20, 000 epe or more. 

The apparent specific heat obtained at each of these frequencies 

(41 

may differ considerably. Kneser' ' has put forth the following ex- 
planation for this phenomenon. As a sound wave of a given frequency 
passes through a gas, a part of the vibrational energy fails to follow 
the acoustic cycle, and the gas departs from its equilibrium condi- 
tion. In fact, a part of the vibrational energy vanishes from the 
expression for the adiabatic elasticity, thus giving an increased 
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sound velocity. Kneser's theory led to the following formula for 
the apparent apccific heat, derived from the velocity of sound data, 

~ M T. X *1. r X. 

C v - Cy f 4 W\ - C,> > <5) 

Cv t- n r Ca 

where C* = true specific heat 

n .» frequency of the sound wave 

S ft 

C* = specific heat of translation and rotation (= for 

a diatomic molecule). 

The values of the apparent specific heats may be substituted in this 
equation to give two equations at each temperature for the determina- 
tion of T and Cy . 

B. IMPACT TUBE METHOD 

In 1942 Kantrowitz^ proposed the impact tube method for 
finding relaxation times and made sufficient measurements to veri- 
fy that the experimental results checked with the theoretical calcu- 

15 ) 

lations that were available at the time. Griffith' later extended 
this work and in 1950 reported on the resulto of his investigations, 
listing the relaxation times determined and expressing confidence 
in the method. 

In the impact tube method, gas at rest is allowed to expand 
adiabatically through a nozzle, and the flow thus produced is brought 
to rest again by compression at the nose of an impact tube (pitot 
tube). The total head defect in the flow, as a measure of the en- 
tropy increase, is directly associated with the heat-capacity lag of 
the gas and can be used to determine the relaxation time. The pro- 
cess is shown schematically in Figure 2. 
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Thc underlying principles involved may be summarized 
briefly as follows. The equations of momentum, continuity, and 
energy for a compressible fluid with possible separate transla- 
tional and vibrational temperatures, together with the standard 
definition for entropy, may be combined to establish the rela- 
tionships, 



4s - [ ^ - Jfl d ( c .ibT,ib') . < 6 > 

d 4- CC,., b T <lb * C P T + (7) 

The notation C y ‘^ refers to the specific heat of vibrational 
motion. JVc have then that »-dE,., b /dT;. >> , and if Cp is the 
specific heat at constant pressure excluding contributions from vi- 
bration, we also have that C|» a . 

( 6 ) 

Using the theory of Landau and Teller' , the variation of 
Cy\V>Tv\b with time may be expressed as 






>TVib) ~ C nb C T ”TV,b) / T 



( 8 ) 



on the basis of considerations involving molecular collisions, the 
relaxation time is found to have the following functional dependence: 

r » p - ' (►'s) \m/T) expfCv's) ] 

The quantities V » s, and m represent, respectively, vibrational 
frequency, molecular radius, and molecular mass. 



The departure from equilibrium between vibrational and 
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translational temperatures was expressed by Kantrowitz as 



€- - C„- Ite (T yIb - T . 



( 10 ) 



In this experiment the changes in T are small compared to T itself, 
and Sq. (6) is reduced to 

jg a fedt T-T v *,b - e'at 

r 






Ul) 



X. 



T T,; b y C.; k T 

T is taken as the average of TT,; b over the whole process. 

The total head defect Ah i3 measured a.3 the pressure dif- 
ference p 0 - , and since this difference is small compared to 

p„ , the differential change in entropy, ds = - Rdp/p, may be closely 
approximated by 



AS * AhR/p„ . 



(12) 



Combining Eqs. (li) and (12) we obtain the relationship, 



Ah/p„ = r . 

J RrC, lb T 



(13) 



If the temperature is eliminated from Eqs. (7), (8), and (10), 
it is found that 



de 

dt 



C^e __ C Vlb <J(qVO 



c ' r 



cL 



dt 



(14) 



For the fast stopping compression process at the nose of the impact 
tube, the second term of this expression becomes small in compari- 
son with the other terms, and the change in e becomes 

Ae - £l . 

C» 2 



( 15 ) 
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Since at the beginning of compression e is zero and the velocity 
is q, , Eq. (15) becomes 



_ C v ; b q, 



C' a 



9 Z 

With q = 0, 2q. (14) is reduced to 

dc 



(16) 



dt 



Sti 

Cpr 



(17) 



For instantaneous compression the result is accordingly, 
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f A* _ e; | 




* J 


2C p RC vi jr j 


l j 



( 18 ) 



At this point it i3 convenient to introduce the following set 
of dimensionless parameters: 

T = Cp<^,T^/ Cpd , e * 2C j 

t' ’ <),t/d , <1 1 l/<|. » 

where d is the diameter of the impact tube. 

If Eqs. (13), (IS), and (19) are now combined, it is found 

that 



(19) 



Ah 



* f t z. 

• - — I e dt . 

;c *' J 



( 20 ) 



In order to determine the value of € in this expression it is neces- 
sary to solve Eq. (14) exactly, the integral being 



€ = exp 





Cfrdt 


h 




rt 4 

Cpdt 


c„ h jCt-) 


j 


't'Cp 




Ydt 


Cp 



dt. , (21) 



which, stated in terms of the dimensionless parameters, leads to 
the expression 



“ exp 



( 22 ) 



-9- 



- 


r 
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r r 1 


r} 
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dt 


exp 





«■ . 
dt' 



The integration of Eq. (22) has been done numerically for 
various values of T # . The integration of Eq. (20) has been car- 
ried out graphically, and from these resulting values a curve of 
Ah/Ah; t vs. T' has been prepared. "With this curve on hand an 
experiment may be conducted using any gas, and the relaxation 
time may be determined from the measured values of T* 0 , p. , 
and pi . The values of C v ;^ , Cp , Cp , and P for the gas can 
be found by the use of well-known statistical methods. 

The procedure then is as follows. ^ 

a) Calculate *T? from the relation ship, T * • 

r V*. 

b) Calculate from the relationship, * |2C P CT 0 -T t M * 

c) Calculate from Eq. (13). 

d) Calculate from the relationship , Ah- p 0 -p 3 . 

e) Use the ratio of dh/Ah lfc to enter into the graph of 
dh/Ahj c vs. T* to find the corresponding value for X 1 . 

f) Calculate the relaxation time, X , from the relationship, 

T - X Cp dy/Cp . 



Yc 



C. EXPERIMENTAL RESULTS 

Table I lists the values of the relaxation times of twenty gases 
determined experimentally by various investigators. The informa- 
tion in this table was obtained, for the most part, from .Reference (5). 
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rn. THEORETICAL. CALCULATIONS 



The accumulation of a large body of experimental data led to 

efforts on the part of investigators to determine vibrational relaxa- 

tion times using theoretical calculations. The method of Landau and 

Teller^, developed in 1936, and refined by Bethe and Teller 4 ^ in 

1940, has received general acceptance. A more detailed method de- 

(8 9 ) 

veloped by Schwarts, Slawsky, and Herzfeld' ’ ’ was published in 
1952. This method of calculation is based on the fact that an equilib- 
rium between the various degrees of freedom can be reached only 
by means of a transfer of energy taking place during a molecular col- 
lision. Thus for a pure gas it can be said that 



T 



(23} 



v/here ^ is the number of collisions a molecule experiences per 
second and is proportional to the density, while P is the probability 
that upon a collision a quantum of vibrational energy will be converted 
to translational energy. For a gas composed of polyatomic molecules 
or mixtures of different types of molecules, there arises the possibility 
of an energy exchange between the different vibrational degrees of 
freedom during a collision as well as exchanges with the translational 
degrees of freedom, and a separate relaxation time must be ascribed 
to each of these energy transfers. 

The theory of Schwarts, Slawsky, and Hersfeld will be used in 
this paper to compute the variation of vibrational excitation probabi- 
lities with temperature of mixtures of O., and N-> » and of HF and H,. 

For a mixture of molecules of type a and b» following the notation of 
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Schwartz et al, we may define the probabilities 



In each of these cases, the molecule, the symbol of which appears 
first in the parenthesis, was initially in the first excited vibrational 
state and finally in the ground or zeroth vibrational state. The sub- 
script designates the behavior of the vibrational quantum number of 
the molecule, the symbol of which appears second in the parenthesis. 
Each of these probabilities has been tabulated and graphed as func- 
tions of temperature for both the O ^ and HF - gaseous mix- 

tures. 

The probability may be obtained from the rela- 

tion 



p 0 ,,(M) = P 0 -,(^ b ) e - Ktv ^/ kT 



( 24 ) 



The others arc determined independently. These probabilities may 
then be used to obtain the relaxation times from the relations 

i 



rS»- ”" P - 






( 25 ) 



and 



L- e“- AT ) • 






( 26 ) 



Here and denote the number of collisions of the designated 

type per second. We may define similarly 'Y (b,b) and Y (b,a). 



a or the complex excitation process we may write 



r 4ib (Ab) " X * X » M «.b P —^' b ) * » 



-W^/kV 



( 27 ) 
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with a similar definition for Here "X^andX^are respec- 

tively the mole fractions of species a and species la in the mixture. 

It should be noted that only under special circumstances can the re- 
laxation times defined above be grouped together to represent an ef- 
fective relaxation time which corresponds to a measured relaxation 
time. 



In general, we may write for the probabilities 




Here ja is the reduced mass of the two colliding molecules, AE. 
represents the amount of energy exchanged between vibration and 
translation, ir 0 and € are the usual constants in the Leonard- Jones 
interaction potential, and Vci) is the customary perturbation 
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integral over harmonic oscillator wave functions. For molecules 
whose value of £ is not too great, it is sufficient to write 
oc*= l7.5/r 0 . In the case of HF, where sufficient virial 
coefficient, or transport data are not available, the collision di- 
ameter, r 0 , was estimated by plotting Vo as a function of the 
electronegativity of the halogen atoms. From Fig. 3, it can be 
seen that the value of Vi, for HF is approximately 1.7A. 

The condition 

| > AE (35) 

will not hold below a certain temperature for each of the probabili- 
ties. Y/e need this inequality to write the expansion 

V f * - V* 1 

It is not necessary that 

although the larger the inequality, the better the convergence. 
Schwartz and Herzfeid, in the appendix of their paper, have given 
a more complex method of computation which applies in the region 
of or below the limit temperatures defined by Eq. (35). 

The following derivation applies for the determination of the 
temperature limit for which the basic equations used here will be 
applicable. The condition required by inequality (35) may be re- 
stated as „ , V/a. 

*o > (^ AE /V/ • 

If we combine this expression with £q. (33) we obtain 



AE (Lit 



f*''? 



V 






♦ - - 



(36) 



» AE 
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/ 4ir*kTAE \ ' y 

\oc^h|A / 



AE 

2.|a \ cO*hp 



m^r * ( 



2 AE 



■)' 






if 



4ir x kAE 



, , - (Ulf • a , 



inequality (37) is reduced to 

(JT) ,/ * * 6y4(jrV / * > 6 , 

■jvhence 

J Vs T Vs “ 6T ,/s - bV4J' /3 >0 . 



(37) 



(38) 



We may obtain from this inequality the lower temperature limit of 
r T* such that 

T iwt - 1-758 e»*/j , 



or 



m 1.24 «t*Vi / AE \ /2 ~ 

-- - TT- hr] • 



(39) 



Eq. (39) has been used to compute the values of the limit tempera- 
tures shown in Table II. 

A compilation of the constants used in making the calcula- 
tions indicated by the basic equations is shown in Table III. The 



constants in the Lennard-Jones interaction potential, and Y m , 



were obtained from data listed in the "Transactions of the AS ME 
with the exception of the value for Y 0 for HF which was estimated 
in the manner already described. The vibrational constant was 
found in each case from "Spectra of Diatomic Molecules" by Herz- 
hers* 11 '. 



,( 10 ) 
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The results of the indicated computations are lie ted in Table 
IV and 3 hown graphically in Figures 4 and 5. The values indicated 
for the Op - mixture agree within a few percent with those listed 
in the paper by Slawsky and Herzfcld. It ic believed that where dif- 
ferences do exist, they may be attributed to a variation in the selec- 
tion of force constants and low-velocity collision diameters inasmuch 
as the published values are not always consistent. The calculated re- 
laxation times also agree reasonably well with values determined 
experimentally, as discussed in the paper of Schwartz et al. No 
data are available for comparison with our results for the HF - Hp 
probabilities. 
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IV. VIBRATIONAL LAG IN THE HYDROGEN-FLUORINE 

ROCKET MOTOR 

It is now possible to present some brief considerations in 
determining whether or not the adiabatic expansion through the 
de Laval nozzle in the - F^ rocket motor may be taken as vi- 
brational near -equilibrium flow. The reactants are taken in such 
proportion that will be in excess, so that the products of com- 
bustion passing through the nozzle will be mostly a mixture of 
HF and 

(12) 

Following a procedure analogous to that of Penner' , in 
his treatment of chemical reaction in nozzle flow, we may take ao 
an estimate fbr the vibrational temperature lag for near -equilibrium 
flow that 

AT — (DT/ot)T . (4U) 

For representative nozzles in 1,000 psia thrust chambers, it has 

been found^*^ that the cooling rate (-DT/Dt) — (Tc-Te^^ 3x10 

°K/sec. The use of very small nozzles will, of course, increase 

the cooling rate above this value; however, for nozzles in larger 

7 o 

thrust engines, (-DT/Dt) will generally be less than 3 x 10~ °K/sec, 
thereby leading to better approximations of near -equilibrium flow 
than we estimate below. 

From Figure 5 it can be seen that the probability P* ..,(*>) 
is approximately as large as some of the other probabilities and 
much larger than others. The smaller probabilities contribute 
little to the effective relaxation time and may be neglected. The 
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probabilities, which are approximately as large as P o e) . 
all contribute to the effective relaxation time. For the purpose of 
simplicity, however, it may be assumed that the maximum possible 
value of the effective relaxation time can be approximated as 

- Ti,b * * (41) 

The actual value of T must, of necessity, be smaller than ae 

the other probabilities, which are approximately as large aa P. 
have not been considered in computing and wouid lower its 

value. If we can show, however, that vibrational near -equilibrium 
flow exists for the calculated upper limit, . the conclusion 

will certainly hold for the real T which is smaller than . 

Since ^,b * which denotes the number of collisions per 
second between HF and molecules, increases markedly with 
pressure, it is apparent that if the conditions at the nozzle exit 
are such that vibrational near -equilibrium flow exists there, then 
it must also exist everywhere in the nozzle. The value of AT 
was therefore computed at the nozzle exit, where it is a maximum. * 

In the present computations, the nozzle exit pressure was 
taken as one atmosphere. The product will generally lie 

between the comparatively narrow limits of 0 . 2.5 and 0. 20, and 
will usually be about 0. 23. The product of ^ and 

n 

will not vary appreciably from 4x10' between the temperatures, 
1400 to 3000° K. Therefore, for exit temperatures within this 
range, £ l x lO sec and, for (-DT/Dt) ^3x10 °K/sec, 

*For a more detailed description of iterative lag calculations dur- 
ing nozzle flow, see Reference (14). 



- 18 - 



* 3* K . 

The estimate that over the length of the nozzle the vibra- 
tional temperature lag cannot exceed 3° K indicates that vibra- 
tional near -equilibrium flow occurs during expansion in a de 
Laval nozzle for representative hydrogen-fluorine rocket mo- 
tors. 
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A. Sonic Velocity Method 

Hydrogen 

Oxygen 

Carbon monoxide 

Chlorine 

Carbon dioxide 

Sulfur dioxide 

Ammonia 

Ethylene 

Methane 

Carbon disulfide 
Nitrous oxide 



B. Impact Tube Method 

i 

Hydrogen 

Nitrogen 

Nitrous oxide 

Carbon dioxide (comm) 

Sulfur dioxide 

Ammonia 

Mythl chloride 

Ethane 

Ethylene 

Propylene 

Propane 

Ethylene oxide 

Butane 

Butadiene 

Methane 



T 1 

°K 


r x 

i 

sees 

i 


293 


i 

i 

0. 018 (rot) j 


294 


1000 


.257 


10 


294 


18 


294 


10. S i 


294 


0. 181 


294 


0. 04 1 


294 


0.233 


383 


0.84 


294 


0.70 


292 


0.92 

! 

i 

! 


287 


\ 

? 

i 

1 

0. 021 (rot) 


289 


0. 002 


291 


1. 12 


291 


2. 02 


293 


0. 37 


290 


0. 12 


295 


Q. 202 


237 


0.003 


289 


0. 207 


294 


0. 006 


294 


0. 007 


303 


1.23 


294 


0. 018 


294 


C. 013 


289 


0.48 



TABLE I. experimental results reported in Reference 5 
for the relaxation times of gases 
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Probability 


Temp ( 


5 k 


°2 * N 2 


HF - H 2 




326° 


1390° 


Po-O <*'.'») 


405° 


2700° 




330° 


2380° 


O C ^*1^) 


403° 


2440° 


Po-*, U,b) 


231° 


592° 




231° 


592° 

i ; 



TAB LE II. Lower temperature limits for which 
the probabilities can be calculated from the theory 
of Schwarts, Slawsky, and Herzfeld. 
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x^robability 


A 


c/k 

®U 


V % to'* 
sac*' 


&E * io‘* 


°2 “ N 2 


n 












1 

P o^o<A<M 


f 

3. 50 


114. 0 


4.75 


3. 14 


P 0-oCb.b) 


3. 69 


93.7 


7. 10 


4. 70 


Po O 


3. 60 


103. 5 


4.75 


3. 14 


P . 


3.60 


103. 5 


7. 10 


4.70 




3. 60 


103. 5 


4.75 

7. 10 


1. 56 


HF - H 0 
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P... (<».«) 
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p.,.^b.b) 
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P. t«.b) 
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8. 22 


P 0 ...<b.«) 


2.32 
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13.2 


8.74 

i 


P.„, (b,«) 


2. 32 

1 


114.0 


12.4 ! .52 

13.2 j 

j ! 



TABLE III. Values of constants used in computations for 
the vibrational excitation probabilities. 
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Variation of probabilities with temperature 
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rig, 1. Schematic drawing of apparatus for determining vibrational 



relaxation time by the sound velocity method 
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TEiif'ERATUltS, HISTORY 




Fig. 2. Schematic drawing of the experiment, ana temperature 

history of tne gas for determining vibrational relaxation time 
by the impact tube method. 
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FIGURE 3- COLLISION DIAMETER AS A FUNCTION 
OF ELECTRONEGATIVITY. 




LOG PROBABILITY 







FlGURE 4- VARIATION OF (a=0 2 ,b = N 2 ) PROBABILITIES WITH 
TEMPERATURE. LIMIT TEMPERATURES CORRESPOND TO 
THE JUNCTION BETWEEN THE SOLID PORTION AND THE 
DOTTED PORTION OF THE CURVES. 
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log Probability 
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FIGURE 5- VARIATION OF (a = HF, b-H^) PROBABILITIES WITH 
TEMPERATURE. LIMIT TEMPERATURES CORRESPOND TO 
THE JUNCTION BETWEEN THE SOLID PORTION AND THE 
DOTTED PORTION OF THE CURVES. 










' 







?33 



P33 



25056 

fenny 

Vibrati nal velaxc tion 
ti >es of gaseous mix- 
tures cf diatoi,iic mole— 
cules and their effect 
on rocket performance. 



25056 

Fenny 

Vibrational relaxation biles 
of gaseous mixtures of diatomic 
molecules and their effect cn 
rocket performance. 



